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Abstract 

Background: Plants respond differently to mechanical wounding and herbivore attack, using distinct pathways for 
defense. The versatile sweet potato sporamin possesses multiple biological functions in response to stress. However, 
the regulation of sporamin gene expression that is activated upon mechanical damage or herbivore attack has not 
been well studied. 

Results: Biochemical analysis revealed that different patterns of Reactive oxygen species (ROS) and antioxidant 
mechanism exist between mechanical wounding (MW) and herbivore attack (HA) in the sweet potato leaf. Using 
LC-ESI-MS (Liquid chromatography electrospray ionization mass spectrometry analysis), only the endogenous JA 
(jasmonic acid) level was found to increase dramatically after MW in a time-dependent manner, whereas both 
endogenous JA and SA (salicylic acid) increase in parallel after HA. Through yeast one-hybrid screening, two 
transcription factors IbNAG (no apical meristem (NAM), Arabidopsis transcription activation factor (ATAF), and 
cup-shaped cotyledon (CUC)) and IbWRKYI were isolated, which interact with the sporamin promoter fragment 
of SWRE (sporamin wounding-responsive element) regulatory sequences. Exogenous application of MeJA (methyl 
jasmonate), SA and DIECA (diethyldithiocarbamic acid, JAs biosynthesis inhibitor) on sweet potato leaves was 
employed, and the results revealed that IbNAG mediated the expression of sporamin through a JA-dependent 
signaling pathway upon MW, whereas both IbNAG and IbWRKYI coordinately regulated sporamin expression 
through JA- and SA-dependent pathways upon HA. Transcriptome analysis identified MYC2/4 and JAZ2/TIFY]0A 
(jasmonate ZIM/tify-domain), the repressor and activator of JA and SA signaling among others, as the genes 
that play an intermediate role in the JA and SA pathways, and these results were further validated by qRT-PCR 
(quantitative real-time polymerase chain reaction). 

Conclusion: This work has improved our understanding of the differential regulatory mechanism of sporamin 
expression. Our study illustrates that sweet potato sporamin expression is differentially induced upon abiotic 
MW and biotic HA that involves IbNAG and IbWRKYI and is dependent on the JA and SA signaling pathways. 
Thus, we established a model to address the plant-wounding response upon physical and biotic damage. 
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Background 

Plants sense and respond to external stimuli using a 
repertoire of mechanisms that regulate gene expres- 
sion for survival in hostile environments. To defend 
against external stimuli, plants have evolved inducible 
defense mechanisms against microbial pathogens and 
herbivores that involve the regulation of gene expres- 
sion for the synthesis of specific proteins and second- 
ary metabolites [1]. Signals that mediate systemic plant 
responses are classified as either slow or fast moving, 
travelling within minutes to several hours, and are typ- 
ically mediated by different hormones/electrical signals 
[2], These fascinating phenomena imply the existence 
of cell-cell communication that transmits the defense 
response over a long distance [3]. Wounding responses 
in plants regulate multiple signaling pathways [4,5]. It 
is well established in tomatoes and Arabidopsis that 
JA is a systemic, long-distance, mobile- signaling mol- 
ecule that functions by transmitting information about 
wounding to distant, non-wounded tissues where a 
defense response is invoked [3,6]. Moreover, jasmo- 
nates, which are synthesized from a-linolenic acid (LA, 
18:3) via the octadecanoid pathway, induce the expres- 
sion of a wide range of defense genes against MW, HA 
and necrotrophic pathogens [2,7]. Mutant and trad- 
itional grafting experiments support that the JA- 
signaling pathway creates long-distance mobile signals 
to activate defense gene expression [8]. In addition, 
there is sufficient evidence that the kinetics of the JA 
and SA signaling pathways aid in choosing the defen- 
sive strategy following stress [9,10]. Cross-talk between 
phytohormones other than JA and SA, ethylene and 
ABA (abscisic acid) has not been extensively studied 
with respect to wounding signal perception (for re- 
views, see [11,12]). It has been well documented that 
signal transduction pathways often overlap and that 
the induced responses to various stimuli are different 
[13]. For instance, insects chewing cause severe damage to 
leaf tissue and release volatile organic compound (VOC), 
which induce direct defense by activating wound signaling 
pathway [5]. Some of the oxylipin compounds (JA) have 
the ability to be a master-switch against herbivores and turn 
on herbivore-related defense genes [14]. Moreover, saliva- 
derived compounds are brought to the wound site during 
HA. These saliva-derived compounds (FAC -fatty acid 
amino acid conjugate) were shown to induce an outburst of 
both JA and herbivore-induced volatile organic compound 
(HI- VOC) [15,16]. How the herbivore feeds on the plant 
tissue, either by sucking or chewing, will determine whether 
the signaling pathway is JA solely, or JA and SA, induced 
[16,17]. The coevolution between plants and herbi- 
vores has generated a costly defensive trait, which is 
greatly dependent on the concentration and timing of 
the hormone (JA, SA, or ethylene) [17,18]. Extensive 



studies defining the global expression profiles upon the 
activation of wounding signaling pathways have been 
conducted using microarrays [19] and proteomic [20]. 
The role of such signaling hormones, integrated in the 
decision-making of stress response, are played by a 
class of Reactive oxygen species (ROS) [21]. An ex- 
ample of ROS as retrograde signaling from chloroplast 
to nucleus having sub-cellular and systemic functional- 
ity has been documented [22,23]. However, these cells 
are equipped with an excellent antioxidant (enzymatic 
or non-enzymatic) defense mechanism to detoxify the 
harmful effects of ROS [24]. Sweet potato (Ipomoea bata- 
tas) with its hexaploid complex genome (2n = 6x = 90) ac- 
counts for the fifth-largest tuber crop worldwide. Unlike 
others, the cultivar Tainong 57 has great potential to defend 
itself against herbivore and/or wounding stress. The versa- 
tility of sweet potato sporamin makes it a very attractive tu- 
berous storage protein because of its strong trypsin 
inhibitory activity and multiple biological functions against 
various stresses [25]. 

In this article, we analyzed the first line of ROS scav- 
enging mechanisms by antioxidant enzymatic activity in 
leaves after MW and HA. We dissect the role phytohor- 
mones play in shaping the interaction between MW and 
HA in the regulation of sporamin expression. We report 
for the first time that the sweet potato leaf, damaged by 
MW and HA, utilizes different phytohormone signals to 
orchestrate the interactive role in the regulation of the spor- 
amin gene. Realizing the potential importance of these 
signals, we studied the regulatory mechanisms that attenu- 
ate two wound-induced transcription factors (IbNACl and 
IbWRKYl) by either JA or SA under the stress of MW and 
HA. Furthermore, we performed a transcriptome analysis 
of wounding stress and surveyed global gene expression to 
identify the relevant genes of the phytohormone network. 
Together, our findings illustrate the precise role of the phy- 
tohormone and signaling pathways that play critical roles in 
regulating sweet potato sporamin expression upon abiotic 
mechanical damage and biotic HA. 

Results 

In vitro antioxidant scavenging enzymatic activities after 
MW or HA 

Initial attempts were made to learn how the different 
ROS levels at the site of damage leads to increasing anti- 
oxidative enzymes and response to stress. Damage caused 
by free radical-induced oxidative stress is repaired by anti- 
oxidant enzymes, which are important as they protect cells. 
We performed two different stresses, MW and HA, to esti- 
mate the total level of antioxidant enzymes. The time 
points were designed to identify the early (0.25 h, 0.5 h) 
and late (6 h) responses of wounding stress in the sweet po- 
tato leaf. Phenolics can directly scavenge molecular ROS to 
eliminate free radicals, linearly correlated with antioxidant 
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capacity. Total phenolic content and 1, l-diphenyl-2-picryl- 
hydrazyl (DPPH) assay was measured for control leaf (cK), 
leaf damaged by MW and leaf damaged by HA. The total 
phenol content from MW leaves increases at 0.25 h (51 ± 
8.02), 0.5 h (45.1 ± 2.3), and 6 h (66.4 ± 11.8) compared to 
that of cK leaves (24.4 ± 5.08) (Figure 1A). The increase in 
phenolic content in response to insect feeding (HA) 
(cK, 24.1 ± 4.42; 0.25 h, 43 ± 1.78; 0.5 h, 41.7 ± 2.2; 6 h, 
15. ±1.05) was likely not significant except for values 
6 h time point (Figure 1A). In addition, the sweet po- 
tato leaves showed a significant increase in DPPH 



antioxidant level in response to MW and HA. Using 
the EC 50 value of ascorbic acid (R 2 = 0.9715), the DPPH 
antioxidant level was 0.90 ± 0.025 ug/mL, 0.36 ± 
0.052 ug/mL and 0.65 ± 0.036 ug/mL at 0.25, 0.5 and 
6 h, respectively, in response to MW (Figure IB). Simi- 
larly, the DPPH antioxidant level was 0.26 ± 0.05 ug/ 
mL, 0.33 ± 0.047 ug/mL and 0.09 ± 0.07 ug/mL at 0.25, 
0.5 and 6 h, respectively, in response to HA. Collect- 
ively, a direct correlation between total phenolic con- 
tent and DPPH free radical scavenging activity was 
shown in response to MW and HA at all the time 




E F 




Figure 1 In vitro antioxidant studies of sweet potato leaves damaged by MW and HA. The third leaves of sweet potato plants were 
sampled at 0.25, 0.5 and 6 h after treatment with MW or HA. (A) The total phenolic content was measured at 725 nm and is expressed in gallic 
acid (GAE) equivalents; (B) The DPPH activity was measured at 517 nm using ascorbic acid as a standard, and the EC 50 was calculated; (C) The 
SOD activity was measured at 560 nm and defined as the amount of enzyme that resulted in 50% NBT inhibition; (D-F) The catalase activity 
(D), POD activity (E) and H 2 0 2 level (F) were measured using the TMB {3,3',5,5'-tetramethylbenzidine) method. Data are means (+SE) of three 
independent experiments with three replicates. {* represents P < 0.05; ** represents P < 0.01; *** represents P < 0.001; 'ns' represents P > 0.05). 
DPPH: 2, 2-diphenyl-1 -picrylhydrazyl; POD: peroxidase; SOD: superoxide dismutase. 
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points tested. However, a significant reduction of the 
total phenolic content and DPPH antioxidant level was 
observed upon HA as compared to MW. 

Antioxidant-catalyzing enzymes (Superoxide dismutase 
(SOD), Catalase (CAT), Peroxidase (POD), and hydrogen 
peroxidase (H 2 0 2 )): The superoxide dismutase (SOD) 
activity increased at all of the time points tested in re- 
sponse to MW and HA, compared to the control (cK) 
(Figure 1C). The mean enzymatic activity of sweet po- 
tato SOD significantly increased at 0.25 h (0.60 U/mL) 
and 6 h (0.59 U/mL) after MW (Figure 1C). However, 
the enzymatic level at 0.5 h (0.30 U/mL) after MW was 
lower than that at 0.25 and 6 h. Similar enzymatic activ- 
ity levels were observed at 0.25 h (0.16 U/mL) and 6 h 
(0.12 U/mL) after HA, but in comparison with MW, HA 
resulted in lower SOD activity (Figure 1C). Catalase 
assay: Catalase (CAT) is an enzyme that reduces oxida- 
tive stress by removing H 2 0 2 . Our results showed that 
in response to MW, the activity of catalase decreased 
significantly over time; 0.5 h (1.42 ± 0.01EA U/mL) and 
6 h (0.75 ± 0.06 EA U/mL) was lower than that at 0.25 h 
(3.2 ±0.45 EA U/mL) (Figure ID). In contrast, upon 
HA, the catalase enzyme activity was significantly re- 
duced at 0.25 h (2±0.19EA U/mL) but increased at 
0.5 h (4.3 ± 0.36EA U/mL) (Figure ID). However, at 6 h 
(0.79 ± 0.06EA U/mL), some activity remained in both 
the MW and HA leaves (Figure ID). Overall, this data 
indicates a significant increase in catalase activity in 
MW and HA leaves compared with control leaves (cK) 
(0.7773 ±EA U/mL). Peroxidase activity assay (POD): 
Hydrogen peroxide production often results in an in- 
crease in peroxidase activity in response to stress. From 
our data, the peroxidase activity did not show any signifi- 
cant differences between the control leaves and the MW 
and HA leaves at the different time points (Figure IE). 
However, a notable increase in the basal level activity was 
observed in the MW leaves. In contrast, the HA leaves 
showed a slight decrease in peroxidase activity at all of the 
time points tested (Figure IE). Hydrogen peroxidase activity 
assay (H 2 0 2 ): Next, we used the TMB (3,3',5,5'-tetra- 
methylbenzidine) assay to determine the hydrogen peroxide 
scavenging activity, which revealed notable changes in this 
activity between control leaves and MW or HA leaves. 
The H 2 0 2 scavenging activity at 0.25 and 0.5 h was 
101.7 ± 1.72 uM and 101.1 ± 2.37 uM upon MW, re- 
spectively, compared to the control (53.3 ± 2.41 uM) 
(1 F). However, the H 2 0 2 scavenging activity gradually 
decreased in both the MW (49.6 ± 2.06 (iM) and HA 
(28.8 ±1.95 uM) leaves at 6 h (1 F). Collectively, our 
results of total phenolic content (Figure 1A), DPPH 
antioxidant level (Figure IB), POD (Figure IE) and 
H 2 0 2 assay (Figure IF) showed similar or identical 
total antioxidant level and activity upon MW and 
HA stresses. In contrast, SOD (Figure 1C) and CAT 



(Figure ID) showed a significantly different pattern of 
total antioxidant level and activity upon MW and HA 
stresses. 

In vivo fluorescence determination of the singlet oxygen 
species ( 1 0 2 ) levels in response to wounding stress 

To assay the singlet oxygen ( 0 2 ) level at an early time 
point is essential because it is the primary signal to in- 
duce other ROS species. Here, we monitored wounded 
leaves (both MW and HA) under the short illumination 
duration (0.25 h) and light intensity. The formation of x 0 2 
was measured by SOSG dye (singlet oxygen sensor green). 
Figure 2A-H shows our confocal laser-scanning microscope 
(CSLM) experiment of leaves from unwounded (cK) and 
MW or HA leaves (0.25 h) that were infiltrated with SOSG 
dye in the dark or light. Substantial SOSG fluorescence 
was only detected in MW or HA leaves treated with 
SOSG in the dark and in the light compared to the cK 
leaf (Figure 2C, D and G, H). The significant gener- 
ation of x 0 2 in sweet potato leaves damaged by MW or 
HA induces hydrogen peroxide and superoxide in the 
presence of light, suggesting the cross-talk role of 1 0 2 
production during ROS metabolism (Figure 2). 

Alteration of endogenous JA and SA levels in sweet 
potato leaves in response to MW or HA 

From the changes of ROS levels, we suspected that 
leaves damaged by MW or HA (Spodoptem littomlis) 
eventually increase endogenous hormones for defense 
signaling. To address this point, a time-course analysis 
(0.25, 0.5 and 6 h) was carried out to quantify the levels 
of endogenous JA and SA by LC-ESI-MS in the third 
leaf of sweet potato plants damaged by MW or HA. In a 
healthy sweet potato leaf, the endogenous levels of JA and 
SA were 79.4 ±4.3 ng/g fresh weight (FW) and 52.3 ± 
1.9 ng/g FW, respectively (Figure 3A). In the leaf damaged 
by MW, the level of endogenous JA began to increase at 
0.25 h and remained significandy higher than its initial level 
at 0.5 h (Figure 3A). This increase in the endogenous JA 
level was 2-fold higher than the control (cK). In contrast, 
the level of endogenous SA began to increase slightly at 
0.5 h and was significantiy reduced at 0.25 h (Figure 3A). 
At 6 h after MW, both endogenous JA and SA were de- 
tected at their basal levels (Figure 3A). In contrast, upon 
HA, the level of endogenous JA positively correlated with 
MW at 0.25 and 0.5 h and at 6 h (2-fold) (Figure 3B). Simi- 
larly, we noticed a significant and stable increase in the SA 
level beginning at 0.5 h (2-fold) and continuing to 6 h (4- 
fold) after HA (Figure 3B). Thus, our results indicate that 
the levels of JA and SA are inversely correlated during 
mechanical damage and HA in sweet potato leaves. The re- 
duction of SA and JA (Figure 3A and B) could be explained 
by the variation between the samples tested. 
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Figure 2 Wound-induced production of ^0 2 in sweet potato leaves. Confocal laser-scanning microscopic (CLSM) images of sweet potato 
leaves. False color image of sweet potato leaf infiltrate with SOSG and exposed to light (150 umol rrf 2 s" 1 PPFD) for 15 min. (A, E) Control leaf 
(cK) without wounding and SOSG treatment that was kept in the dark for 1 2 h. (B, F) Control leaf infiltrated with SOSG and kept in dark for 
15 min. (C, G) Wounded leaf (MW, HA) infiltrated with SOSG and exposed to dark for 15 min. (D, H) Wounded leaf (MW, HA) infiltrated with SOSG 
and exposed to light for 15mins. (I, J) Relative quantification of image inflorescence by Image J. The circle denotes the wounded area. The black 
arrow indicates the production of ] 0 2 . The SOSG fluorescence was measured using an excitation wavelength of 525 nm. SOSG: singlet oxygen 
sensor green; W: wounded; cK: non-wounded. 



NAC and WRKY were the DNA-binding proteins that 
recognize SWRE element of sporamin promoter 

A yeast one-hybrid screening system was performed to 
isolate a DNA-binding protein that binds to three tandem 
repeats of the SWRE fragment {3xSWRE (38-bp), -1103- 
ACATTTCTCGTAAATACGTACAATATCCTTGTCTT 
TCC -1061). A cDNA library generated from wounded 
sweet potato leaves was expressed as a translational fu- 
sion with a GAL4 activation domain (AD/cDNA library) 
in a yeast reporter strain which carries an integrated 
His3 allele, with the three tandem copies of the SWRE 
sequence (3xSWRE) employed as a probe (Figure 4A). 
At the first screening using the library, which contains 
6 x 10 cfu, 94 colonies were selected from the selective 
medium and were re-streaked on a more stringent 



selective medium. Finally, two cDNAs, named IbAMCl 
(Accession: GQ280387) and IbWRKYl (Accession: 
GQ280386) respectively, emerged within the SWRE frag- 
ment sequence. In order to confirm their binding ability 
with SWRE fragments, these two clones were transformed 
into a yeast cell containing the pHIS2-3xSWRE reporter 
vector or pHIS2 as a negative control again. As shown in 
Figure 4B, co-transformation of pGADT7-Rec2-VoNAC\ 
with pHIS2-3xSWRE or pHIS2 can both grow in the SD/- 
LeuATrp medium, as well as pGADT7-Rec2-IbWRI<Yl 
with pHIS2-SWRE or pHIS2. However in the SD/-His/- 
Leu/-Trp plus 100 mM 3- AT selective medium, only 
co-transformation of pGADT7-Rec2-lbNACl or pGADT7- 
Rec2-lbWRKYl with pHIS2-3xSWRE can grow but not 
pHIS2 reporter vectors only (Figure 4). 
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Figure 3 LC-ESI-MS quantification of the endogenous 
phytohormones JA and SA in sweet potato leaves. The healthy 
third leaf of a sweet potato plant without treatment was used as a 
control (cK). The JA and SA levels were quantified at 0.25, 0.5 and 
6 h after the third leaf was damaged by: (A) MW; or (B) HA (second 
instar larvae of Spodoptera littoralis). The level of JA and SA was 
expressed as ng JA/SA per g FW. (* represents P < 0.05; ** represents 
P < 0.01; *** represents P < 0.001; 'ns' non-significant represents P > 
0.05). Values are the means (±SE) of three replicates (n = 6). 



The regulation of sporamin is facilitated by transcription 
factors in response to wounding stress and hormonal 
treatment 

To address whether the systemic expression of spora- 
min, IbNACl and IbWRKYl gene were activated by 
endogenous hormones JA and SA in response to wound- 
ing, we executed an experiment in which similar time 
points were analyzed. QRT-PCR was performed on 
sweet potato leaves treated subjected to MW and treated 
with the exogenous application of MeJA, SA or DIECA. 
We reproducibly demonstrated that the expression of 
sporamin transcript is induced in leaves damaged by 
MW over the time course (Figure 5A). Similarly, the ex- 
pression of IbAMCl was strongly induced at 0.25 and 
0.5 h and at 6 h (Figure 5B) after leaves damaged by 
MW. 

On the other hand, the accumulation of the IbWRKYl 
transcript did not show any induction at all of the time 



points tested (Figure 5C). Thus, sporamin and IbAMCl 
were induced upon MW stress, suggesting the associ- 
ation of IbAMCl as a positive regulator and IbWRKYl as 
a negative regulator of sporamin gene expression upon 
MW stress. 

Next, we examined the expression level of the spora- 
min, IbNACl and IbWRKYl gene after the application 
of MeJA. Using qRT-PCR, we confirmed that the spora- 
min and IbAMCl transcripts induced at 0.25 and 0.5 h 
after the exogenous application of MeJA, and their level 
noticeably declined from 0.5 to 6 h (Figure 6A and B). 
Furthermore, IbWRKYl was not detected at the early 
time points (0.25 and 0.5 h) and showed a remarkably 
higher expression at 6 h after MeJA treatment (Figure 6C). 
In plants, DIECA is a potent inhibitor of JA- biosynthesis 
in the octadecanoid pathway [26]. As expected, exogen- 
ous application of DIECA was only effective against the 
expression of sporamin transcript at the early time 
points 0.25 and 0.5 tested, and had no effect at 6 h after 
treatment (Figure 6D). However, exogenous application 
of DIECA does not alter the IbAMCl gene expression 
(Figure 6E). In contrast, IbWRKYl transcript had no ex- 
pression at early time points, similar to the without- 
DIECA treatment (Figure 5C), but induced to express 
only at 6 h after DIECA treatment (Figure 6F). Together, 
this can be easily explained by the fact that blocking of 
JA biosynthesis will automatically lower the signal 
transduction in the JA-signaling pathway; antagonistic- 
ally, SA levels should be induced and SA-responsive 
genes should be up-regulated. 

Furthermore, upon treatment with SA, our qRT-PCR 
results shows that the low level or basal expression of 
sporamin was significant at all the time points tested 
(Figure 6G). On the other hand, strong activation of 
IbAMCl gene expression was stable post-exogenous 
application of SA treatment at all the time points 
(Figure 6H). Thus, we hypothesize IbNACl is a posi- 
tive convergence regulator, activated upstream of sporamin 
that is influenced by both JA- and SA-mediated signaling in 
a time-dependent manner. Moreover, the accumulation of 
IbWRKYl did not show any expression at early time points 
(0.25 and 0.5), but induced to express at 6 h after SA treat- 
ment (Figure 61). Therefore, IbWRKYl has a partial role 
in both JA and SA signaling at later time points in the 
wound- signaling pathway. 

Effects of HA on the expression of sporamin and the 
transcription factors lb/MCI and foWRKYl in sweet 
potato leaves 

To demonstrate that leaves damaged by HA have differ- 
ent responses to those with MW, second instar larvae 
(starved for 2 days) were placed on the underside of 
every third leaf of sweet potato plants and the leaves 
were observed for biotic damage. Herbivore {Spodoptera 
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Figure 4 Isolation of IbAMCI and \bWRKY'\ proteins interacting with SWRE DNA by yeast one-hybrid screening. Two cDNA clones 
(IbNACl and IbWRKYI) were isolated by Y1H screening, and re-transformed into the yeast strain Y187 in the presence of plasmid containing three 
tandem repeats of the SWRE sequence (3xSWRE-pHIS2) (A) or without 3xSWRE3 (p/-//52). Growth of yeast transformants is shown on both + His 
(A) and -His/100 mM 3-AT medium (B). Positive control and negative control are employed to standardize selection condition by using yeast 
strain Y187 co-transformed with pGAD-Rec2-53/p53HIS2 or pGAD-Rec2-53/pHIS2 (Clontech). 



littoralis)-eaten leaves were collected at three different 
time points (0.25, 0.5 and 6 h), and the mass of the 
leaves consumed by the larvae was moderately reduced 
at 6 h after feeding began (Figure 7A). We further ob- 
served the behavior and survival of the herbivores, and 
our results showed that the body color of the larvae 
changed from brown to brownish-green and the larvae 
moved from the upper side to the lower side of the 
leaves whenever perturbed. Furthermore, the growth of 
the larvae was measured every 2 days for a period of 
7 days. A few of the larvae died at various stages, leaving 
6-8 individuals that survived to adults (Figure 7B). Their 
weight moderately increased, beginning at day 4 and 
continuing to increase until day 7 (Figure 7C). The lar- 
vae produced wet fecula from the dorsal anus from day 
1 to day 3 (data not shown). 

Next, we examined the expression of spommin, IbAZACl 
and IbWRKYI in response to HA by qRT-PCR. We ob- 
served a similar pattern of spommin and IbAMCl gene 
expression between MW (Figure 5A and B) and HA 
(Figure 8 A and B). However, we did find a significant 
accumulation of IbWRKYI transcripts in leaves dam- 
aged by HA (Figure 8C). In the HA leaf, the accumula- 
tion of IbWRKYI began at 0.25 h, peaked at 0.5 h and 
was maintained until 6 h (Figure 8C). This differs from 
leaves with MW (Figure 5C), in which the transcript 
levels of IbWRKYI transcription factor was inversely 
expressed. Collectively, we concluded that endogenous 
JA plays a major role in the activation of spommin 
gene expression upon two different stresses (MW or 
HA) mediated by IbNACl. However, the differentiated 
responses of lb WRKY through S. littomlis-speciiic oral 



secretion created a cross-talk interaction of JA and SA 
signaling pathways. 

Solexa sequencing: global expression profiling of the 

sweet potato response to MW 

Digital gene expression library sequencing 

In this study, we conducted Solexa sequencing from sweet 
potato leaf (third leaf) with MW (15mins) and control 
without wounding (cK). A total 41,237,622 million raw 
reads were produced by using the llumina Solexa Genome 
Analyzer (GA) II sequencing platform, with total nucleo- 
tides 2, 784,289,408 (10 MB). The mean read length, contig 
length and GC % were 75.00, 519 bp and 44.74% for con- 
trol and 76.15, 493 bp and 44.70 for wounding library. 
Distribution of total clean reads and its coverage are shown 
(Additional file 1: Figure SI). After quality trimming, the 
removal of primers, adaptor, poly (A) sequences, and low 
quality reads, we obtained a total of 20,777, 585 and 
16,098,635 reads from MW and cK libraries, respectively 
(Additional file 2: Table SI). A total of 36,876,220 reads 
were obtained from the wounding and control libraries 
(Additional file 3: Table S2). A total of 41,806 differentially 
expressed transcripts were identified from control and 
wounding libraries. We have analyzed the most differen- 
tially regulated transcripts with a log2ratio >1 or < -1 as a 
threshold by means of significant value (P < 0.001) as well 
as false discovery rates (FDR < 0.01), representing 1,070 dif- 
ferentially expressed (684 up- and 386 down-regulated) 
transcripts (Figure 9) (Additional file 4). The major 
GO (Gene Ontology) terms corresponding to the dif- 
ferentially expressed transcripts are shown in Figure 10. 
A large number of transcripts were associated with 
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Figure 5 QRT-PCR analysis of the sporamin, lb/MCI and IbWR/CVI 
transcript after different wounding treatments. Relative expression of 
sporamin, \bNAC1 and \bWRKY1 expression were measured by qRT-PCR 
in untreated (cK) and mechanically wounded (MW) sweet potato leaves 
at the indicated time points. Ubiquitin (UBQ) gene was used as 
an internal control. (* represents P < 0.05; ** represents P < 0.01 ; 
*** represents P < 0.001; 'ns' not significant represents P > 0.05). 
The results are representative of three replicates. 



'response to stimulus', 'signaling', 'hormone response', 
'metabolic' and 'cellular process'. Among the differen- 
tially expressed genes (DEGs), we observed that a large 
number of transcripts were associated with responses 
to stimulus such as allene oxide synthase (AOS), alter- 
native oxidase (AOX), sensitive to proton rhizotoxicity 
(STOP1), SEN1, salt tolerance zinc finger (STZ), non- 
race specific disease resistance (NDR1) and nucleosome 
assemble protein- related protein NRP1, including spora- 
min (SPOR), ipomoein (IPO), etc. Furthermore, we no- 
ticed that the several transcription factors, such as 
Arabidopsis transcription activator factor (ATAF2), 
WRKY, Arabidopsis response regulator (ARR), ethylene 
responsive factor (ERF), MYC etc., overlap in the sub- 
category of 'response to stimulus'. 

Enriched pathway analysis of DEGs 

In this study, the DEGs were mapped to a total of 89 Kyoto 
Encyclopedia of Genes and Genomes (KEGG) pathways 
with Q value of <0.05 (Figure 11). The 15 most enriched 
pathways were involved in 'phenylalanine metabolism', 
'starch and sucrose metabolism', 'nitrogen metabolism', 
'plant-hormone signal transduction', 'amino sugar nu- 
cleotide sugar metabolism', 'plant-pathogen interaction', 
'phenylpropanoid biosynthesis', 'alanine, aspartate and 
glutamate metabolism' and 'biosynthesis of secondary 
metabolites'. Altogether, these results shed light on the 
wounding response in the sweet potato leaf, which is 
overwhelmed by dynamic metabolism as well as actively 
adapting to the hostile environment (Figure 11). 

Gene enrichment analysis 

Next, GO term enrichment was analyzed using BiNGO 
software (http://bingoware.sourceforge.net/). We next 
screened and compared the DEGs against the Arabi- 
dopsis genome database. Of the 1,070 differentially 
expressed transcripts, 809 (497 up- and 312 down- 
regulated) transcripts were matched to the Arabidopsis 
genome with an E-value 10' 10 . Of these 809 differen- 
tially expressed transcripts, 590 genes were found to 
have putative biological functions (Additional file 3: 
Table S2). These 590 genes were assigned to various 
functional role categories based on their significant 
matches to proteins that were already assigned a functional 
role. Their expression patterns were separated into six 
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Figure 6 Expression profiles for sporamin, lbW>4C1 and IbW/WVI in sweet potato leaves. The third leaves of sweet potato plants were 
sampled at 0.25, 0.5 and 6 h after treatment with MeJA (0.5 uM) (A-C), DIECA (50 mM) (D-F), or SA (2 mM) (G-l) independently. Total RNA 
was extracted at the indicated time points from untreated (cK) and treated leaf samples. (A-l) Relative expression of sporamin, IbMACI and 
\bWRKY] transcripts were measured by qRT-PCR. Asterisks indicate statistically significant differences (Dunnett's multiple comparison test, 
Represents P < 0.05; ** represents P < 0.01; *** represents P < 0.001; 'ns' not significant represents P > 0.05). The results are representative of 
three replicates. 



major clusters based on their gene enrichment values. 
Transcripts from each cluster are presented in the 
Additional file 5: Figures S2 and Additional file 6: Figure S3. 

Quantitative real-time PCR (qRT-PCR) confirming accuracy 
of Solexa sequencing 

The expression of all 15 DEGs was consistent with the pre- 
dictions based on the Illumina sequencing results. Among 
them, nine genes (alpha-glucan phosphorylase 2; UDP-D- 
glucuronate 4 epimerase 1; Ethylene-responsive transcription 
factor 3; Tetraspanin8, OBP3-responsive proteinl; Zinc 



finger CCCH domain containing protein 29; carbonic 
anhydrase 1; late embryogenesis abundant hydroxyl 
proline-rich glycoprotein; and copper transport family 
protein) were up-regulated and six genes (expansin-like B2 
precursor; wound-responsive protein; histone Hl.l; short- 
chain dehydrogenase/reductase (SDR) family protein; and 
EFE and GAMMA-VPE) were down- regulated (Additional 
file 7: Table S3). Actin2 was chosen as the reference gene 
for normalization data. In most cases, the expression values 
were a little higher than those obtained from the Illumina 
sequencing (Figure 12). 



Rajendran et al. BMC Plant Biology 2014, 14:1 12 
http://www.biomedcentral.com/1471-2229/14/112 



Page 10 of 21 



cK 0.25h 0.5h 6h 



B 




Spodoptcra littoral is 



Insect feeding bioassay 



£ 0.4 

OS 

«'¥ 0-3 

5 * 



'a 2 
o > 0.1 




infestation 



I) 



Spodoptcra littoralis 



| 3 0.15 
B JS 

u o.io 
> o 
< S 

C-.05 



0 2 5 

Days post inoculation 



Days post inoculation 



Figure 7 Insect bioassay. (A) The mass of leaves consumed by herbivores {Spodoptera littoralis) at the indicated time points; (B-C) The mean mass (±SE) 
of second instar larvae of Spodoptera littoralis reared on sweet potato leaves (n = 10) and the plant growth biomass were recorded after 7 days (B) and 
from 2-7 days (C) of feeding (D). The results are representative of three replicates. (* represents growth impairment of larvae). 



Early wounding-responsive genes that act in the first 
line of the wounding signaling cascade as assessed by 
qRT-PCR 

Based on our computational analysis, we have selected 
genes related to the plant-pathogen interaction path- 
way and plant-hormone signal transduction, and we 
discovered 10 differentially expressed, early wounding- 
responsive genes which may act in the first line of the 
wounding signaling cascade, including cyclic nucleotide- 
gated ion channel (CNGC), calmodulin (CaM), JA-induced 
WRKY, MYC2, MYC4, JAZ1, JAZ4, JAZ6, non-expressor of 
pathogenesis related (NPR1), and TGA. We compared the 
time-dependent changes in the expression of these genes 
with qRT-PCR (Figure 13). The results showed that JAZ1/ 
TIFY 10b, JAZ2/TIFY 10a, MYC2 and MYC4 were sig- 
nificantly up-regulated at early time points and these 
up-regulated genes were associated with wounding or the 
JA- mediated signaling pathway. Therefore, the induced ex- 
pression of JAZ1/TIFY 10b, JAZ2/TIFY 10a, MYC2 and 
MYC4 confirms that the first line of the signaling cascade 
requires JA signaling in the control of the defense response. 



On the other hand, the down-regulation of CNGC 
and CaM in response to wounding suggests their roles 
as negative regulators of defense gene expression in 
the wounding response. It is possible that Ca + ICaM 
control downstream defense genes in the wounding 
signaling pathway. In addition, we examined NPR1 and 
TGA, which were down-regulated in response to wound- 
ing, which suggests the NPR1- and SA-dependent stimulus 
pathway (Figure 14). 

Discussion 

To minimize the damage caused by either MW or HA, 
plants produce antioxidant and ROS scavenging en- 
zymes [27,28]. Different patterns of ROS and antioxidant 
mechanisms exist between MW and HA in the sweet 
potato leaf (Figure 1). This study supports the notion 
that two different stresses (MW and HA) generate ROS 
in the plant cells. During this process, production of 1 0 2 
initiates lipid peroxidation to disrupt PSII as a result of 
plasma membrane NADPH oxidase (Figure 2) [29,30]. 
The enhanced activity of SOD in response to MW 
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Figure 8 HA and gene expression profiles for sporamin, \bNAC\ 
and \bWRKV\ in sweet potato leaves. (A-C) Relative expression of 
sporamin, IbMACI and IbWRKY] transcripts were measured by qRT- 
PCR analysis in untreated leaves (cK) and leaves damaged by HA. 
The results are representative of three replicates. Asterisks indicate 
statistically significant differences. (Dunnett's multiple comparison 
test* represents P < 0.05; ** represents P < 0.01 ; *** represents P < 0.001 ; 
'ns' not significant represents P > 0.05). The results are representative of 
three replicates. 

» J 

constitutes the primary defense that efficiently protects 
membrane lipids against oxidative damage [31]. In paral- 
lel, MW reduced CAT activity that can be explained by 
its interaction with SA, and induced the role in plant- 
pathogen defense mechanism [28]. However, the increase 
in CAT activity upon HA may prevent excess H2O2 pro- 
duction, which is responsible for membrane lipid peroxi- 
dation [32]. Increasing evidence of differential activation 
of these antioxidant enzymes was solely dependent upon 
growth and development [33]. 

Sweet potato has differentially and chronologically 
expressed antioxidant enzymes when under MW and 
HA stress. The antioxidant system protection against 
wounding stress depends on growth stages. SOD acts 
as a primary defense signal for MW at all time points 
tested. Conversely, CAT and POD followed a two-step 
activation upon MW and HA. At an early time point, 
CAT was deactivated upon MW stress. In this instance, 
accumulation of hydrogen peroxide acts as a signal to 
induce octadecanoid signaling transduction pathway. 
In contrast, CAT is activated after prolonged stress 
caused by HA, and leads to a decomposition of excess 
H 2 0 2 at an early time point. At the later time point, 
both CAT and POD are deactivated. Accumulated 
H2O2 acts as a signal to induce octadecanoid and other 
defense signal transduction pathways (Figures 1 and 2) 
[33-36]. Furthermore, our study presented evidence 
that a specific level of JA and SA production induces 
the expression of defense-related genes that are in- 
volved in the antioxidant and detoxification processes 
(Figure 3) [28,37,38]. 

Systemin and JA-mediated signal in sweet potato leaf 
after MW and HA 

We investigated the critical role of phytohormones in 
orchestrating sporamin gene expression in response to 
MW and HA. The key observation underlying the exper- 
iments described is that endogenous JA level was en- 
hanced after MW on the leaf to mediate the signal 
transduction and trigger the activation of the octadeca- 
noid pathway (JA pathway) during early wounding time 
points (Figure 3) [39,40]. However, upon HA, the en- 
hanced endogenous SA level could be due to larval oral 
secretion (violicitin), which interacts with endogenous 
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Figure 1 1 KEGG pathways of DEGs responding to wounding stress. 



phytohormones to elicit the SA burst and attenuates 
JA (Figure 3) [40,41]. In agreement with previous stud- 
ies of sporamin induction after MW or HA (Figure 4) 
[42,43], our report presented evidence that the sporamin 
is wound-inducible and activated through a systemin- 
mediated JA-signaling pathway [44]. Although JA is neces- 
sary for sporamin expression in both MW and HA, while 
SA was shown to be antagonistic to sporamin expression 
in sweet potato (Figure 5). These differences in cross-talk 
were consistent with previous reports of MW and HA in 
leaves, as they are solely dependent on timing, concentra- 
tion of endogenous signal, and the regulation of defense 
gene expression. [16,17,39,45-47]. 



Wounding signal transduction is mediated by two 
transcription factors to regulate sporamin gene 
expression differentially 

Many wound-responsive promoters in plant species 
have been functionally characterized, and several puta- 
tive wound-responsive cis-elements have been identi- 
fied such as G-box [48], GCC box [49], AG-motif [50,51], 
Jasmonate/elicitor responsive element [52,53], 13-bp/L-box 
[54], W-box, [53] and sporamin [42]. In the present work, 
we describe the isolation and characterization of two novel 
transcription factors that bind within the 38-bp -1103- 
ACATTTCTCGTAAATACGTACAATATCCTTGTCTTT 
CC -1061- of the SWRE region of the sporamin promoter 
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Figure 12 Validation of transcriptome data by qRT-PCR. Genes 
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stress (0.25 h) condition. Quantitative real-time PCR analysis was 

performed for 15 genes, and the values (least squares means) were 

fitted to a linear regression model. Both the x- and y-axes are shown 

in a log 2 scale. Real-time PCR data was obtained from three independent 
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of the sweet potato (Figure 4). Furthermore, we identified 
signal cross-talk between JA and SA signaling pathways 
that were mediated by IbAZACl and IbWRKYl differen- 
tially (Figures 5 and 6). It appears that minor variations in 
the core sequences of sporamin promoter impart respon- 
siveness to these different stimuli. Upon MW, IbNACl, 
acting as an early transcriptional factor, functions in the JA 
signal pathway to regulate sporamin gene induction in the 
sweet potato leaf (Figures 5B and 6B) [55]. However, 
IbNACl may partially function in the SA pathway, which 
harmonizes the action of other downstream defense genes 
(e.g., PR genes) (Figure 6H) [56]. Similarly, induction of 
lbNA CI expression upon HA explains its role as a conver- 
gence transcriptional activator, functioning in both JA and 
SA signaling pathways to regulate sporamin gene expres- 
sion in the sweet potato leaf (Figure 8B) [55,56]. In our 
study, we were puzzled to find that IbWRKY is not a posi- 
tive regulator for sporamin gene expression in MW or SA 
signaling (Figures 5C, 6G and I). However, the MeJA/SA 
-mediated inhibition of VoWRKY expression at an early 
time point suggested that endogenous JA level might play 
an important role in controlling its expression (Figure 6C 
and I). Moreover, the expression of IbWRKYl after 6 h 
upon exogenous MeJA or SA treatment (Figure 6C and I) 
was a consequence of re-methylation of JA for a positive 



feedback effect on SA biosynthesis [42,57-60]. This can be 
tested by quantifying the individual metabolites of JA and 
SA biosynthesis, and remains to be elucidated. 

Interestingly, the strongly up-regulated expression of 
IbWRKYl following HA (Figure 8C) suggested, in this 
instance, at least two steps of regulation were imposed: 
firstly, an early response of IbWRKY, possibly under the 
control of S. littoralis-specific oral secretion that elicits 
SA burst. This holds true because the feeding by S. lit- 
toralis caused more of a marked increase of endogenous 
SA than the endogenous JA concentration (Figure 3). 
Secondly, the impact of leaves continuously damaged by 
S. littorais, different from 'single' wounding events, accu- 
mulates VoWRKYl expression possibly by FACs, through 
which JA level was maintained [60,61]. In light of this 
consideration, VoWRKYl played a direct role in the 
defense against herbivore and, indirectly, maintained JA 
level for sporamin- activated defense. 

To make a comparison for the cross-talk differences 
between the two transcription factors IbNAC and 
IbWRKY in JA and SA signaling, we inferred from the 
publication and evidence by Li et al. (2004) that the 
inactivated WRKY70 upon JA signaling to the promoter re- 
gion of JA- responsive gene was likely controlled by cyto- 
solic NPRlprotein. Needless to say, IbWRKY could be 
analogous to WRKY70 as they reflect a similar cellular SA 
and JA balance. Therefore, the fate of IbWRKY regulates 
the sporamin gene expression via integrating SA and JA- 
signaling events, depending upon the stress. In support of 
this assumption, IbWRKY and WRKY70, an uncommon 
member of the WRKY transcription factor, are pivotal inte- 
grators of the SA and JA signals to regulate plant defense 
responses. 

In contrast, our IbNAC is strictly responsive to both 
wounding stresses (MW or HA) involved in both JA- and 
SA- signaling pathways, which suggests it is a positive tran- 
scriptional activator for regulating sporamin gene expres- 
sion and partially other downstream genes (e.g. PR), Upon 
prolonged HA feeding, NPR1 is not analogous to IbNAC 
but acts as an amplification loop for IbWRKY and PR gene 
expression through the SA-mediated independent pathway 
to defend against additional entry of pathogen signal [57]. 
In conclusion, our results support the need for future work 
on investigating the mutants of IbAMCl and IbWRKYl and 
their respective roles. 

Global gene expression profiles in sweet potato response 
to MW stress 

Based on the previous literature, MW stress often over- 
lapped with multiple stress responsive genes and shared be- 
tween different plant species [13,62-64]. The sweet potato 
is a non-model organism that has not been extensively 
studied. Among the contig annotated, we found the major- 
ity of DEGs involved in the 'response to stimulus' and 
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Figure 13 Temporal expression patterns of the DEGs in sweet potato leaves damaged by MW. Relative gene expression measured by 
qRT-PCR at 0.25, 0.5 and 6 h, comparing untreated leaves (cK) and leaves damaged by MW. (A-E) Up-regulated genes: JA-induced WRKY, MYC2, 
MYC4, JAZ1-TIFY-10A, and JAZ2-TIFY-10B. (F-J) Down-regulated genes: CNGC, CaM, NPR1, TGA, and JAZ6-TIFY-11B. Error bars indicate the standard 
deviations for the means of triplicate samples, and the Act2 gene was used for data normalization. Asterisks indicate statistically significant 
differences. (Dunnett's multiple comparison test,* represents P < 0.05; ** represents P < 0.01; *** represents P < 0.001; 'ns' not significant 
represents P > 0.05). The results are representative of three replicates. 



'response to stress' were related to pathogen, wounding, 
and abiotic and biotic stresses (Additional file 5: Figure S2 
and Additional file 6: Figure S3). Based on the KEGG path- 
way (Figure 11), we found that the majority of metabolic 
pathways genes were enriched and are common for sweet 
potato, suggesting they are involved in cell wall biosyn- 
thesis, inter-conversion of nucleotide sugars, biosynthesis of 
lignin and other cell-wall-bound phenolic compounds at 
the site of leaf damage caused by MW [65,66]. On this 
basis, we integrated the results from this study and other 
published studies and formed a working model for the acti- 
vation of spommin gene expression and its regulation in re- 
sponse to external stimuli (Figure 14). MW or HA activates 
JA production, which in turn independendy activates the 
expression of IbAMC and spommin through a JA-mediated 
signaling pathway. Once the JA signal reaches its threshold 
level, genes encoding JA biosynthetic enzymes, which 
are tightly regulated by JA-responsive transcription 
factors (TFs), are up- or down-regulated. The activa- 
tion of IbNACl binds to the SWRE and induces transcrip- 
tion of the spommin gene during early response of 
wounding stress. The negative cross-talk role of VoWRKYl 



specifically functions in SA- signaling and has an indirect 
role in JA- signaling pathway for the persisting expression 
of spommin. 

In contrast, HA results in an increased biosynthesis of 
both SA- and JA- signaling and leads to the subsequent 
activation of parallel signaling pathways. In this context, 
both IbNACl and IbWRKYl coordinate with each other 
and bind to the SWRE to activate the expression of 
the spommin gene for defense. Under prolonged HA, 
IbWRKYl possibly interacts with FACs, which is ne- 
cessary for maintaining the early action of de novo JA 
biosynthesis in the leaves. Therefore, IbWRKYl is in- 
volved in both pathways upon HA, and requires the 
initial signal strength (either extremities of high or 
low) for the mutual antagonistic interaction of the SA- 
and JA-mediated signaling to regulate spommin gene 
expression for defense. 

Furthermore, the selected genes involved in the first 
line of wounding signaling cascade support the partial 
cross-talk between JA- and SA-mediated signal pathways 
in the regulation of spommin gene expression. For cor- 
relation, the spommin promoter region also has a G- 
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Figure 14 Proposed model for sporamin gene regulation in sweet potato leaves. Two mechanisms induced by wounding were clarified in 
this work: activated by abiotic MW or biotic HA. Our findings illustrate that both of them have an impact on different ROS levels and antioxidant 
mechanisms, which indicates signal components are different. The results presented two wound-induced signals: MW or HA can change the 
endogenous level of JA or SA to regulate sporamin. MW was shown to be JA-mediated, signaling to regulate the transcription factor IbNACI 
for the activation of sporamin gene expression. On the other hand, HA was shown to be both JA- and SA-mediated, signaling to regulate 
transcription factors both IbNACI and IbWRKYI for the activation of sporamin against herbivores. This antagonistic difference was further 
confirmed through the analysis of genes identified in the transcriptome and biological network in the context of MW stress. The selected 
genes involved in the first line of signaling cascade supports further the partial cross-talk between JA- and SA-mediated signal pathways 
in regulation of sporamin gene expression. 



box-like element and GC-core-like element, indicating 
that there is room to directly target MYC/JAZ interac- 
tions. Similarly, the binding site of W-box for IbWRKYI 
was controlled by calmodulin through other TFs (e.g., 
WRKY/TGA) of downstream defense genes. Thus, the 
wounding-induced expression of sweet potato sporamin 
differs between MW and HA due to the a's-elements of 
the sporamin promoter. However, this is regulated by 
IbNACI and IbWRKYI and depends on the JA and SA 
pathways (Figure 14) [67]. 

Conclusion 

Two mechanisms induced by wounding were clarified in 
this work: differentially activated by abiotic MW or bi- 
otic HA. Our findings illustrate that both of them have 
an impact on different ROS levels and antioxidant mech- 
anisms, which indicates that signal components are dif- 
ferent. The results showed that the two wound-induced 
signals, MW or HA, can change the endogenous level of 
JA or SA to regulate sporamin. MW was shown to 
be JA-mediated, signaling to regulate the transcription 
factor IbNAC for the activation of sporamin gene ex- 
pression. On the other hand, HA was shown to be both 
JA- and SA-mediated, signaling to regulate transcription 
factors of both IbNACI and IbWRKYI for the activation 
of sporamin against herbivores. This antagonistic differ- 
ence was further confirmed through the analysis of 



genes identified in the transcriptome in the context of 
MW stress. The selected genes involved in the first line 
of signaling cascade supports further the partial cross- 
talk between JA- and SA-mediated signal pathways in 
regulation of sporamin gene expression. The knowledge 
gained from this study holds great potential for spora- 
min research within the wounding signaling pathway in 
the near future. 

Methods 

Plant material and growth conditions 

Sweet potato (Ipomoea batatas cv. Tainong 57) plants 
were grown in pots at 28°C under a 16-h light (2000 
lux)/8-h dark (16 L/8D) cycle. Sweet potato with at least 
five fully expanded leaves per pot (n = 3) were utilized 
for each experiment, and each data set was from three 
biologically independent replicates. 

Antioxidant enzyme extraction and activity assay 

The third sweet potato leaf was collected at various time 
points (15 min, 30 min and 6 h) after MW and HA and 
frozen immediately in liquid N 2 (as described above). 
For enzyme assay, 0.3 g of leaf was ground with 3 mL of 
ice-cold 50 mM sodium phosphate buffer (pH 7.5). En- 
zyme assays were then performed immediately. Total 
phenolic and DPPH radical scavenging assays were per- 
formed spectophotometrically [68,69]. The extraction 
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buffer for estimating total phenolic content contained 
0.5 mL of sodium carbonate (5% w/v) and the solution of 
gallic acid (0-0.1 mg/mL) in ethanol (95%). For DPPH, the 
extraction buffer contained the solution of ascorbic acid 
(5-25 ug), 0.1% DPPH solution [69]. The extraction buffer 
for SOD (EC 1.15.1.1) contained 100 uL of 50 mM potas- 
sium phosphate (pH 7.8), 0.1 mM Na 2 EDTA, 14 mM- me- 
thionine, 75 uM -NBT and 2 uM - riboflavin [70]. The 
extraction buffer for CAT (EC 1.11.1.6) contained 0.6 mL 
potassium phosphate buffer (pH7.4), 0.5 mL hydrogen per- 
oxide (H 2 0 2 ) (molar extension coefficient: 44 M' 1 cm' 1 ) 
[71]. For POD (EC1.11.1.7) and H202 assay, the extraction 
buffer contained 0.05 mL of guaicol solution (20 mM) and 
0.03 mL of 12.3 mM H2O2. Hydrogen peroxidase activity 
assay was performed using the TMB method. The enzyme 
extract buffer contained 20 mM TMB dissolved in 95% 
ethanol, 600 uL of 0.1 M acetate/citrate buffer (pH 4.5), 
100 uL of 100 uM of H202, 100 uL of 0.5 M H 2 S0 4 and 
100 uL [34,72]. 

Fluorescence detection of singlet oxygen species 
generation in response to wounding 

Fully expanded, intact third leaves were wounded by 
MW or HA at 0.25 h. SOSG (Invitrogen) was dissolved 
in methanol at a final concentration of 50 urn in 
50 mM sodium phosphate buffer (pH 7.4), according 
to the manufacturer's instructions. Unwounded (cK) 
and wounded leaf was slowly infiltrated with 50 um 
SOSG either in the dark or light (150 umol m' 2 s' 1 
PPFD) for 30 min. The leaves were dried in the dark 
and then photographed using CLSM with excitation at 
475 nm and emission at 525 nm [73]. 

Phytohormone analysis of JA and SA by LC-ESI-MS 

Unwounded (cK), MW and HA sweet potato 3rd leaf 
were sampled at 0.25, 0.5 and 6 h and frozen immedi- 
ately at -80°C. The samples (100 mg) were immediately 
ground to a fine powder in liquid N 2 and then exposed 
to extraction buffer (100 ul of 100% methanol). The 
samples were centrifuged with 1 mm and 2 mm (1:1) 
zirconium oxide beads by vortexing at 16,000 x g for 
10 min. After centrifugation at 16,000 x g for 10 min, 
the supernatant was collected for LC-ESI-MS analyses 
based on the procedure [74]. 

Construction and screening the one-hybrid library and 
screening interacting proteins 

A total 950 ug of RNA from sweet potato 3rd leaf 0.5 h 
and 1 h after wounding was used to extract poly (A) + 
RNA by Oligotex mRAN Midi Kit (QIAGEN). Briefly, 
aliquots of poly (A) + RNA were reverse-transcribed separ- 
ately by using either random or oligo (dT) primers, and the 
cDNA products were combined and inserted in the shuttle 
vector, pGADT7-Rec2, containing GAL4activation domain 



(pGAL4AD). Yeast one-hybrid reporter (Matchmaker™ 
One-Hybrid Library Construction & Screening Kit) 
construct was prepared by cloning tandems containing 
three repeats of cis-elements of sporamin promoter 
{3xSWRE (38-bp)-1103- ACATTTCTCGTAAATACG- 
TACAATATCCTTGTCTTTCC -1061) was used as a 
bait sequence, and inserted in reporter plasmids pHIS2 
to select the DNA-binding domain encoded in the 
sporamin cDNA library. The pHIS2 reporter vector of 
3xSWRE, cDNA and pGADT7-Rec2 vectors were co- 
transformed to the yeast strain Y187 lacking His-, 
Trp-, and Leu-. The transformed yeast cells were spread on 
SD/His/-Leu/-Trp in the presence of 40 mM 3-AT. Trans- 
formation efficiency was carried out by LiCl- polyethylene 
glycol, and transformants grown on SD/-Leu, SD/-Trp were 
tested for ji-galactosidase activity. Plasmids from putative 
positive clones of 300-400 colonies appeared on SD/-His/- 
Leu/-Trp +40 mM 3-AT plates were isolated after three 
days of incubation. The 252 colonies were analyzed by 
colony PCR to amplify the cDNA insert. Ninety-four 
amplified cDNA inserts were sequenced, among them 
two cDNA with putative domains as conserved, namely 
NAC and WRKY. To eliminate false positive colonies, 
these plasmids were separately introduced to the yeast 
cells containing either 3xSWRE or pGADTl -Rec2- 
IbNAC/IbWRKY grown in the SD/-His/-Leu/-Trp + 
100 mM 3-AT selection medium. 

MW, application of hormone and hormone inhibitor and 
herbivore feeding experiment 

The wounding of the leaf was carried out as follows. The 
third leaf counted from the top of the apical meristem 
was wounded by squeezing with forceps (one time) and 
then collected after 0.25, 0.5 and 6 h. Plants with fully 
expanded six to eight leaves per pot were sprayed with 
the following chemicals: MeJA, SA and diethyldithiocar- 
bamic acid (DIECA). The chemical and concentrations 
were used as follows: MeJA (0.05 mM) dissolved in 75% 
ethanol, SA (2 mM) dissolved in water and DIECA 
(1 mM) dissolved in water. A third leaf was collected at 
0.25, 0.5 and 6 h after treatment. For control, tissues 
were sprayed with water. All the chemicals were pur- 
chased from Sigma. For insect feeding, second instar lar- 
vae of Spodoptera littoralis were starved for 2 days and 
subsequently allowed to feed on the sweet potato leaves. 
Approximately 7 to 10 second instar larvae per pot were 
placed on the upper side of the third intact leaf and were 
observed until they began feeding [75]. The larvae that 
failed to feed on the leaves within 30 min were dis- 
carded. The leaves eaten by the larvae were collected at 
different time points, 0.25, 0.5 and 6 h, and stored 
at -80°C until use. Larvae were weighed on an analytical 
balance. Feeding was allowed for 7 days, and the masses 
of the leaves consumed by the larvae and masses of the 
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larvae were recorded for every other day for 7 days. Each 
trial had a total of 10 second instar larvae (n = 30) and 
was repeated three times. Mortality and longevity of the 
larvae were recorded daily. On day 7, all the live larvae 
were weighed and mortality was recorded. 

QRT-PCR 

Total RNA was prepared from sweet potato leaf tissues 
using the pine tree method. [76] The quality of the RNA 
was assessed using Tecan Nanodrop analyzer (Tecan 
group Ltd, USA, San Jose). The A260/A280 ratios of the 
RNA samples, which were dissolved in 10 mM Tris 
(pH 7.6), ranged from 1.9 to 2.1. Total RNA (5 ug) was 
used to synthesize cDNA following the instructions of 
the First Strand cDNA Synthesis Kits (Fermentas). 
Completely synthesized cDNA was utilized for RT-PCR 
using specific primers. The primer sets used are listed as 
follows: SPOR F'-ATTAGTTTGGACGGGCG', R-'TTAAC 
GCCGCCGTTAAGG'; IbNACl (Accession:GQ280387) 
(F'CGGCCGGGGATACAAATTTGTAAGCTT, R-'GAAT 
CGGAATCCCGGCGGCATCTC); VoWRKYl (Accession: 
GQ280386) (F-'GGAATTGAGATGGCTGTAGAGTTGT 
TG',R-'GAAAACCCACTAACTTGCGGCTGAGA) and 
Ubiquitin (UBQ) F-('CAACAACACTGCCGTCGC, R-TC 
TTCTCCTTCATCATCGCCAC'). 

Illumina cDNA library preparation and sequencing 

Total RNA was extracted from each leaf sample using a 
previously described Cetyl trimethylammonium (CTAB) 
procedure (see above). The A260/A280 ratios of the 
RNA samples, which were dissolved in 10 mM Tris 
(pH 7.6), ranged from 1.9 to 2.1. For cDNA synthesis 
and Solexa sequencing, 20 |ig total RNA was used at a 
concentration of >400 ng/ul. Poly(A) mRNA was first 
purified using oligo(dT) beads then the mRNA was frag- 
mented into small pieces (100-400 bp) using divalent 
cations at an elevated temperature. The cleaved RNA 
fragments were then primed with random hexamers and 
used to synthesize first-strand and second-strand cDNA. 
The cDNA fragments that were >200 bp in size were 
selected by agarose gel electrophoresis and enriched by 
PCR amplification. Finally, two single-end cDNA libraries 
were constructed for sequencing on an Illumina (GA) IIx 
platform. 

Pre-processing and quality assurance of sequencing reads 

The reads delivered by the Beijing Genomics Institute 
(BGI) were processed using the following procedure to ac- 
quire qualified reads. First, the reads containing adapter se- 
quences or including more than 10% uncertain nucleotides 
(Ns) were excluded. Next, a read was further excluded if 
50% of its bases were low quality (quality score <5, Illumina 
GA Pipeline vl.3). 



De novo assembly 

The qualified reads of wound-stress and control samples 
were respectively assembled into distinct contigs using 
SOAP denovo software (http://soap.genomics.org.cn), which 
applies the de Bruijn graph data structure to construct con- 
tigs (assembled sequences). All of the parameter settings 
followed the standard procedures suggested by BGI. Next, 
the two sets of contigs were integrated into one 'unigene 
set'. The sequences of this 'unigene set' were used as the 
reference sequences in the subsequent analysis of gene 
quantification [77]. 

Gene annotation of the assembled sequences 

The assembled sequences were annotated using a previ- 
ously described procedure involving sequential BLAST 
searches (Altschul, Madden, Schaffer, Zhang, Zhang, Miller 
& Lipman, 1997), which are designed to find the most de- 
scriptive annotation for each sequence suggested by stand- 
ard BGI analysis. The assembled sequences were compared 
with protein sequences collected from GenBank's non- 
redundant database [http://www.ncbi.nlm.nih.gov/RefSeq/] 
using the BLASTx algorithm. In the BLASTx analysis, the 
e-value cutoff was set as le~ 10 . 

Quantification and differential expression analysis of genes 

All filtered reads were mapped to the reference se- 
quences and then the accumulated read counts were 
normalized as the number of mapped reads in a given 
contig per kilobase of contig length per million mapped 
reads (RPKM) [78]. To identify the genes stimulated by 
wounding stress, the expression values of each gene in 
the two samples were estimated. The statistical signifi- 
cance of the differential expression of each gene was 
determined using the method described [79], and the re- 
sults of the statistical tests were corrected for multiple 
testing with the Benjamini-Hochberg FDR. Sequences 
were deemed to be significantly differentially expressed 
if the adjusted P value was <0.001 and there was at least 
a two-fold change (>1 or < -1 in the log2 ratio value) be- 
tween the two samples. All of the procedures followed 
the standard BGI analysis. 

Pathway and GO enrichment analysis 

The KEGG, the major public pathway-related database, was 
employed to identify the metabolic or signal transduction 
pathways containing the DEGs from the whole genome 
background using a hypergeometric test [80]. The same cal- 
culation was applied for GO enrichment analysis and both 
analyses were incorporated into the standard BGI analysis. 

Validation of DEGs by semi-quantitative and quantitative 
real time-PCR (qRT-PCR) analysis 

Total RNA was isolated from control and wounded (15- 
min) samples of fully expanded third leaf tissue. First, 
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we performed semi-quantitative RT-PCR analysis for 15 
genes and submitted them for re-sequencing analysis to 
confirm the reliability of the Illumina sequencing. Subse- 
quently, 15 genes were subjected to validation by qRT- 
PCR. The temporal expression pattern of the remaining 
10 genes was subject to study at various time points 
as described above. All primer pairs are listed in the 
Additional file 8: Table S4. Quantitative RT-PCR was 
performed with an ABI 7500 Real-Time PCR system in a 
final reaction volume of 20 ul containing 5 ul diluted 
cDNA in water, 10 ul 2x SYBR Green Real-Time PCR 
Master Mix (A) and 10 uM each of forward and reverse 
primers. The thermal cycling conditions were as follows: 
40 cycles of 95°C for 3 sec for denaturation and 60°C for 
45 sec for annealing. QRT-PCR was performed for three 
biological replicates. Actin2 (Act2) was used to normalize 
gene expression. The relative expression levels of the 
genes were calculated using the 2 AACT method [81], 
which determines the difference in the CT between the 
control Act2 product and the target gene product. Rela- 
tive fold change for each gene was compared to control, 
which was set to one. 

Statistical analysis 

All statistical analyses were performed using GraphPad 
Prism 5.0 (GraphPad Software, La Jolla, CA, USA). Bio- 
chemical analysis, the quantification of phytohormones, 
was performed by two-way analysis of variance (ANOVA) 
(Bonferroni post- test). One-way ANOVA followed by 
Dunnett's Multiple Comparison Test at 95% confi- 
dence intervals) was performed for gene expression. 
Bar represents mean (SE±). 
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